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Abstract 

We propose the multiple reflection expansion as a tool for the calcu- 
lation of heat kernel coefficients. As an example, we give the coefficients 
for a sphere as a finite sum over reflections, obtaining as a byproduct 
a relation between the coefficients for Dirichlet and Neumann boundary 
conditions. Further, we calculate the heat kernel coefficients for the most 
general matching conditions on the surface of a sphere, including those 
cases corresponding to the presence of delta and delta prime background 
potentials. In the latter case, the multiple reflection expansion is shown to 
be non-convergent. 



1 Introduction 

Heat kernel coefficients play an important role in many areas of theoretical 
physics. They govern the short- distance behavior of the propagator and the 
small-time asymptotics of the Schrodinger equation. In quantum field theory, 



* e-mail: Michael . Bordag@itp . uni-leipzig. de 
T e-mail : vassil@itp . uni-leipzig. de 
te-mail:falomir@obclix. fisica. unlp.edu. ar 
§e-mail:mariel@obelix. fisica. unlp.edu. ar 



heat kernel coefficients define the one-loop counter-terms and quantum anomalies, 
as well as the large mass expansion of the effective action. It is clear, therefore, 
that it is important to have an effective method of calculation of these coefficients. 

To the best of our knowledge, heat kernel methods were first applied to quan- 
tum physics by Fock in 1937 |L|; then, they were re- introduced by Schwinger in 
the 50's (see 0). Due to De Witt ||, these methods became standard in quan- 
tum field theory. The De Witt iteration procedure proved to work quite well on 
manifolds without boundaries and (after certain improvements) allowed to cal- 
culate many terms in the asymptotic expansion of the heat kernel [Q, f| H, Q. 
On manifolds with boundaries, the methods based on functorial properties of the 
heat kernel || ||, [HJ appeared to be more appropriate. These methods allowed to 
calculate some higher terms of the heat kernel expansion, e.g, for local boundary 



conditions [11] and for a transmittal problem [12]. The functorial methods being 



the most general and the most powerful ones, they still have their limitations. 
They work particularly well for most general operators in a certain category. 
However, the number of independent invariants which can enter a heat kernel 
coefficient grows very fast with the order of the asymptotic expansion, so that 
combinatorics becomes unmanageable. Alternatively, the general Seeley calcu- 
lus, which is applicable for general boundaries, may be used. But this method 
becomes unwieldy beyond low orders. 

The methods mentioned above are analytical, i.e., they produce local formulae 
for the heat kernel coefficients in terms of the relevant geometric invariants. An 
alternative to such methods are special case calculations (see e.g. O, ITJ, 



16 , [P7|l and references therein). In this case, the complexity of the calculations is 
almost independent from the order of the asymptotic expansion, but the method 
can only be applied to those problems where a sufficiently high symmetry allows 
for the separation of variables. 

Another alternative is provided by iterative, resp., recursive methods. Well 
known is the DeWitt iteration method. Less known are re-formulations in terms 
of integral equations. For example, in the Lipmann-Schwinger equation for 
the scattering problem was used to determine the asymptotic expansion of the 
Jost function entering the regularized ground state energy. While the known 
iterative methods work well for sufficiently smooth background fields, an effec- 
tive method working for singular background fields or for boundary conditions is 
missing. 

It is an aim of the present paper to suggest the multiple reflection expansion as 
such method. In fact, it is based on an integral equation whose kernel is located 
on the boundary. The iteration of this equation gives rise to the multiple reflec- 
tion expansion. The important point is that only a finite number of reflections 
contribute to a given heat kernel coefficient. 

The use of the multiple reflection expansion in connection with vacuum energy 
is not new. In |19| it was employed to investigate the asymptotic density of 
eigenvalues which, in the modern language, is equivalent to the calculation of heat 
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kernel coefficients. In [20] the possibility of using the multiple reflection expansion 
was mentioned, but found to be too complicated for a general boundary. As far 
as we know, the method has never been used as a tool for the calculation of the 
heat kernel coefficients. In this paper, we demonstrate its effectiveness with a 
recalculation of the heat kernel coefficients on a sphere. As a nice byproduct, we 
re-obtain a representation of the coefficients as finite sums where the difference 
between Dirichlet and Neumann boundary conditions resides in certain signs only. 

In general, the multiple reflection expansion can be viewed as some kind of 
perturbative expansion. For instance, for imaginary frequencies it provides a well 
convergent series for the propagator. It should be mentioned that, in certain cases 
(for instance, with Dirichlet boundary conditions ), this is so despite the absence 
of a small expansion parameter. So, the question is whether this convergent 
behaviour is a general feature. The answer is no, and we provide a counterexample 
by considering the most general matching conditions on the surface of a sphere. 
They are described by a four-parameter family and correspond, for instance, to 
the presence of a delta function or its derivative on the surface. It turns out that 
there is no expansion in powers of the parameter in front of the derivative of the 
delta function but, instead, a nice expansion in the inverse of this parameter, 
which cannot be at all obtained by a multiple reflection expansion. 

We would like to note that a Green function with a ^'-function potential has 
been considered before, for example, by path integral methods in ||21|| . There, 
it was noticed that a perturbative expansion similar to that for a 5-function 
potential yields some unsolvable relations, a fact that is not surprising in view of 
the non-analyticity found by us. 

The paper is organized as follows: In Sec. 2 we collect the necessary formu- 
las on spectral functions and their relations to the heat kernel coefficients. We 
write down the multiple reflection expansion in terms of integral equations for the 
propagator as well as for the heat kernel. In Sec. 3 we use the multiple reflection 
expansion in order to re-obtain the heat kernel coefficient for the classical bound- 
ary conditions on the sphere. In Sec. 4 we consider the most general matching 
conditions on a sphere, and calculate the corresponding heat kernel coefficients. 
Sec. 5 contains a discussion of our results and the conclusions. Some useful 
formulas are banned into Appendix 1, while Appendix 2 contains the study of 
matching conditions in higher dimensional spaces. 



2 Spectral functions 

In this section, we define the spectral functions to be used in the rest of the present 
paper, and give a short introduction to perturbative methods, supplemented with 
some examples of application. 

Let us consider the Laplace operator on a domain Q 6 M. D . Let $ n (x) be 
its eigenfunctions, fulfilling Dirichlet or Neumann (or, more generally, Robin) 
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boundary conditions on 5 = dfl, X n being the corresponding eigenvalues : 

- A$ n (x) = A„$ n (f). (1) 

(In Sec. § we will consider the more complicated case of matching conditions on 
a surface in IR D ). We consider three local spectral functions. The first one is the 
resolvent (propagator) D w (x,y) (at imaginary frequency), obeying the equation 

{u 2 -A)D u3 (x,y) = 6(x-y). (2) 

It can be represented as 

D „ (f|S) = E *^w. (3) 

UJ -f- A n 

n 

The second spectral function is the zeta function, given by 

((x,y;s) = J2K s ®n(mn(y) (4) 

n 

and, finally, the third is the heat kernel 

K(x,y\t) = ^$„(f)$ n (y)e- iA ". (5) 

n 

These functions are related by means of 

POO 

T(s)((x,y;s) = / dt f^K&flt) 
Jo 



2 f°° 

dwu^D^x.y). (6) 

,1 - 8) Jo 



In addition, we consider the corresponding global quantities, which appear 
as integrals over Q of the local ones in the coincidence limit. Because of the 
distributional character of the heat kernel coefficients a n (x,x), it is useful to 
introduce a test function f(x) into this integration. So, let 

C[/](s)= f dx f{x)((x,x;s) resp. K[f](t) = / dx f{x)K(x, x \ t) (7) 
Jn Jo, 

be the global zeta function (resp. heat kernel). In many cases (as, e.g., for mani- 
folds with boundaries and local boundary conditions) the latter has an asymptotic 
expansion as t I 

K m)~ jAur, E a n{f\t n - (8) 

1 ^ n=0,i,l,... 



We should warn the reader that the existence of the expansion (|j) cannot be 
taken for granted. For example, in the case of some pseudo-differential operators 
or non-local boundary conditions, In t terms appear. 

If the expansion (H) exists, one can take / = 1 to define the global heat kernel 
coefficients: 

a n = a n [l) (9) 

If, apart from an appropriate behavior at large t, the heat kernel has a power 
law asymptotics at small t, the zeta function C[/]( s ) is a meromorphic function 
of s. From Eqs. (|7|) and the coefficients a n can be represented by the 
corresponding residua: 

a n [f]= Res (47r)^ 2 r( S )C[/]( S ) (n = 0,1,1,...). (10) 
s=f-n 1 

Furthermore, we remind the reader that, in general, the coefficients consist of a 
bulk and a surface contribution 



fln[/] = / dxf(x)b n (x) + / dfj l (a)f(a)c n (a), (11) 
Jn Jan 

where we have used a G dQ as a notation for a point on the boundary, as opposite 
to the genuine notation x G Vt for a point in the bulk. 

Now, we integrate Eq. (|6|) over the domain Q and insert the result into Eq. 
fllPP. We thus arrive at 

a n [f] = ^s \ } / duj ^- 2s / dx f{x) D w (x,x) (12) 
s = f -n F ( 1 ~ s > Jo Jn 

as the basic equation for calculating the coefficients out of the propagator. 

Having briefly reviewed some basic definitions and well-known facts, we pro- 
ceed now to a brief presentation of perturbative methods. 

The perturbative expansion for the resolvent is constructed in the following 
way: Let D°(x, y) be a zeroth order resolvent. Usually, D° is taken to be the free 
propagator in a flat space without boundaries. Consider the Dyson equation 

D(x , y) = D°(x, y) + J dzD°(x, z)LD(z, y) , (13) 

where the integration goes over a sub-manifold E C £1, and L is some operator 
associated with the perturbation (see examples below). Equation ( |13~D has the 
formal solution: 

D(x, y) = D°(x, y) + J2j dzi - J dz n D°(x } z^LD ^) . . . LD°(z n , y) . (14) 



n=l s 
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In [52] it was shown that the heat kernel has a similar representation, 



oo n 2 

K(x,y;t) = K {x,y;t) +^2(-l) n J ds n J ds n - X ...j ds x J dz n . . . J dz x 

n=1 dM dM 

xK (x,z n ;t- s n )LK ( ) . . .Li^ (5*i,y;si), (15) 

where K Q is a suitable chosen "zeroth order" heat kernel. 

In order to clarify these definitions we consider, in what follows, some exam- 
ples. Let D° = (-V)- 1 be the propagator for a second order differential operator 
T>. Let £ = Q and let L be multiplication by a potential V(x). This is the 
standard situation of the DeWitt expansion with a smooth background potential 
written in form of an iterated integral equation. Then the equations ([13]) - ( p~5| ) 
follow from the formal expansion of D = (— V + V)' 1 and K = exp(— 1(— T> + V)) 
respectively, so that D(x, y) and K(x, y\t) are the propagator and the heat kernel 
of the operator (—T> + V). If the potential V is smooth and falls off sufficiently 
fast at infinity, all integrals in (fL5|) exist. From dimensional considerations, it 
is clear that the highest power of V which may contribute to the heat kernel 
coefficient a n is V n . Therefore, only the first n terms of the expansion (HD must 
be taken into account. 

In our next example, let the operator L be again the multiplication by a 
potential V, but now, let S be a subsurface of co-dimension 1 in In [p2 ], it was 



shown that all terms in the expansion (|15|) exist and give a power-law asymptotics 
of the heat kernel. Later, this expansion was used in actual calculations of the 



heat kernel coefficients [E3~ 



The Dyson equation is useful also for rather general perturbations of boundary 
conditions as, e.g., for the case where more derivative terms are added to the 



usual Neumann one (see [24]). In this case, however, dimensional arguments do 
not work and an infinite number of terms contribute to any given heat kernel 
coefficient a n . 

These examples demonstrate that the "common sense" arguments work rather 
well. If there is a parameter e in the theory such that there is a smooth limit 
e — > of the heat kernel coefficients (such as V — > above), then the formal 
expansions fllil) , fll5|) in that parameter usually give a good approximation for 
the spectral functions. If such parameter is of positive mass dimension, only a 
finite number of terms contribute to each a n . 

This is, however, not the end of the story. In the next section we will see 
that one can construct a perturbative expansion, the so called multiple reflection 
expansion, even when no parameter or limiting procedure exist. Moreover, also in 
this case, only a finite number of terms contribute to each heat kernel coefficient 
- a result which is hard to predict on the basis of "common sense" arguments. 



1 This problem is a particular case of a more general transmittal problem (see Sec. 
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3 Multiple reflection expansion applied to the 
heat kernel coefficients for the sphere 

This section contains a short overview of one particular perturbative method, 
which is particularly well-suited for the treatment of boundary problems, i.e., 
the multiple reflection expansion for Dirichlet and Robin boundary conditions. 
Balian and Bloch applied this expansion, in the boundary value problem 



context, to calculate the density of eigenvalues, which is related to the heat 
kernel by a simple integral transformation. In their work |[25|| , they pointed out 



that the divergent part of the Casimir energy is given by a few first reflection 
contributions. This fact, however, has not been fully appreciated. Therefore, 
we find it useful to repeat some basic facts, translating them to a more modern 
language, and to supply the reader with a simple example. In doing so, we omit 
many details which can be found in the original literature |19[ ^] . 



The multiple reflection expansion is based on simple formulas known from 
electrostatics: Having in mind the application to Dirichlet boundary conditions, 
let /i(a) be the density of a double layer (dipole layer) on a surface S. The 
corresponding potential is 



$(£) = / d 2 a x A^(f-ai) d Sl //(ai), 
Js 



(16) 



where d g 1 is the normal derivative, restricted to the surface S and acting to the 
left. Explicitly written, it reads A u (x — a) da = n(y)V yA^x — y)\ g=s , where n 
is the normal vector. The measure on 5* is d 2 a\ = du\du2^fg where (wi,^) are 
the coordinates of a point a{u\,U2) on S and gij = J^J^f- is the metric. In Eq. 
(|ToT), the propagator A^ix — y) is the free one, i.e., without boundary conditions. 
In three dimensions it is simply the Yukawa potential 

A UJ (x-y) = - (r=\x-y\). (17) 

477 r 

The potential <3>(af) is discontinuous for x approaching the surface 5* (x — > a) and 
the equation 

limg^ s ^(x) = I d 2 a x A w (a - a x ) ^ + -/J,(a) (18) 

Js ^ 

holds. The additional contribution (last term) appears due to the fact that limit 
and integration do not commute. 

In a similar fashion, having in mind the application to Neumann boundary 
conditions, the potential x(5f) of a charged surface with charge density pia) 

x{x) = I d 2 ai A w (x - ai) p{di) (19) 
Js 
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has a discontinuous derivative, 



lim nV x x{x) = [ d 2 a x ~B \ A w (a-ax) p(ai) - -p(a). (20) 
x-^d J s Z 

In general, the multiple reflection expansion for the resolvent reads 

Du,{x, y) = A w (f - y)+K / d 2 a 1 A UJ (x - a x ) dsi A UJ (a 1 - y) (21) 

Js 

+k 2 / c/ 2 ai / d 2 a 2 A LO {x - a x ) d 3l A aJ (ai - a 2 ) <9a 2 A w (a 2 - y) 
Js Js 

+ . . . 

with the notation d= d + d . For k = 1 this propagator obeys Dirichlet and, 
for k — — 1, Neumann boundary conditions. The validity of this expression can 
be verified by noting that it fulfills the differential equation for x ^ S. Moreover, 
boundary conditions can be checked using (^) and (^), whereby the additional 
contributions give rise to cancellations between successive orders of reflections. 
The expansion (|2l| ) is called multiple reflection expansion because it can be inter- 
preted as a motion described by the free propagator from x to being reflected 
(however under any angle due to the integration over a x ), further moving to a 2 
and so forth. More details can be found in |]26| and related papersQ. 



A simple example for the multiple reflection expansion appears if the surface 
S is a sphere. In this case, the expansion becomes an algebraic one. It can be ob- 
tained from fl2T|) by turning to spherical coordinates. It is, however, easier to use 
the known expression for the exact propagator with given boundary conditions, 

DUx,x') =^Y l , m {0,<p)YZ m (8>,t/)D l (rrf, (22) 



l,m 



with 



D,(r, r') = -= [I^ur^K^ury) - I v {ur)I v {ur')K DjR ) (23) 
yrr 



and v = I + Here, I u (x) and K„(x) are the modified Bessel functions, and we 



2 

have introduced the notations 



for Dirichlet boundary conditions , and 



^ (R u KJujR)) 
Kr = a g J' (25) 



8R 

& (JW„(wi2)) 



2 It must be stressed that despite its simple form, the derivation of equation ( ^l| ) contains 
several subtle points, which are explained in the Appendix of pq|. 
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for Robin boundary conditions, where the solutions of Eq. ([!]) have to fulfill 
d_ / r «+5 n ( r ) j = 0. For u = — |, these reduce to Neumann boundary condi- 

tions on the 2-dimensional sphere. 

The multiple reflection expansion appears in the following way |26| . Represent 

Kd ~ I v (uR)KUuR) ' (26) 

and use the Wronskian I , v (x)K u (x) — I v (x)K' v {x) = 1/x to rewrite the denomi- 
nator in ( PSD as 

I v {luR)K'„{uR) = ^L(\- ujR^ {I v {uR)K v {uR))\ . (27) 
Next, expand this denominator so that one obtains for Kjj the representation 



d x 1 



duR 



(I u (ujR)K v (ujR)) . (2 



OO s 

K D = -2ojRK u {uR)K' v {ojR) ^ luoR 

k=0 ^ 

In a similar way, one obtains 

K R = -2ujRK v {ujR) (K(ujR) + -^-K(uR)) (29) 

V ujR / 

x £(-i) M( ^ + :j^M)^M) 

J,— n V \ 



This formal expansion has been shown [gfj to be equivalent to the multiple re- 
flection expansion (|2l|), where the number of reflections is k + 1. 

In view of Eq. (12), we perform the integration over the surface of the sphere 
and define 



D u (t)= / d/i(3)D w (f,x) = 5](2Z + l)A(r,r) 
Jan i=0 



(30) 



so that the coefficients a n , Eq. (pj), turn out to be given by 



Res 4^ 7 / drr 2 V(2/+l) / du u^'DAr, r). (31) 

3 „ r(l - s) Jo ^ Jo 



-R 



S =2~ n 



The procedure to calculate the coefficients from this representation is as fol- 
lows: First, we remark that the poles in s result from large, both uj and I, in the 
Bessel functions. The poles corresponding to boundary contributions (the c n in 
Eq. (0)) appear, in addition, from the upper limit of the integration over r. So, 
we use the uniform asymptotic expansion of the Bessel functions (it is given in 
the Appendix), together with the multiple reflection expansion ( p8|) or ( p9|) , and 
insert them into Eq. fl3~T|). 



Let us start with the first contribution in the rhs. of Eq. fl2"3"|). It corresponds 
to the free space propagator and, thus, it doesn't know about the boundary. 
Consequently, it gives the volume contribution, which is ao = ^-R 3 . 

In order to calculate the higher coefficients, we consider the second term in 
the rhs. of Eq. (|23|) . According to the sum in the rhs. of Eq. fl28"|) resp. Eq. ( f29[) 
we represent the coefficients as a sum over reflections, 



2n 
fc=0 



(32) 



Using the uniform asymptotic expansion of the Bessel functions these coefficients 
can be calculated (for details, see Appendix 1). As a result, for Dirichlet boundary 
conditions, the first 

&ri S £1X6! 



k= 







1 


2 


3 


4 


5 




n = 


l 

2 


-2tt 3 / 2 












-2tt 3 / 2 


n = 


1 


2tt 


h 










1- 


n = 


3 
2 








-1^3/2 
6" 










n = 


2 









-¥ 






315 


n = 


5 
2 













48 


12" 


-^7T 3 / 2 
120 


n = 


3 





40 71 
3003 


-lis 71 " 




— lis 71 " 




9009 



The coefficients a n are proportional to R 3 ~ 2n . In this table we have taken 
R = 1. When replaced in Eq. (|3^), they sum up to the known values (shown 
in the last column) as can be checked, for example, by comparing them with the 
results in appendix B of [15|]). It is interesting to note that the corresponding 
heat kernel coefficients for Robin boundary conditions with u = can be obtained 



(k) 

through the replacement a n — > {— iJ'On 

As a last example in this section, we give some reflection contributions to the 
heat kernel coefficient a 2 for Robin boundary conditions. They read (R = 1) 



.(0) 



a? 

a? 
J 3 ) 



0. 

105 

7T 

105 

7T 

315 



'-3 + 12u - 28u 2 



(5 - 42m + 140m 2 - 280m 3 ), 
(35 - 270m + 756u 2 - 840m 3 



The sum (|32|) gives the known result a 2 = — 18m + 60m 2 — 120m 3 ) 

particular, Neumann boundary conditions on a sphere follow by choosing u = 



In 

_i 

2 - 
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4 Singular potentials on a spherical shell 



In this section we will study the heat kernel expansion for the free Laplacian in 
IR D , acting on the space of functions obeying on a D — 1-dimensional sphere, 
S D ~ 1 , certain matching conditions relating the values of the functions and their 
first derivatives on different sides of the sphere. If one assumes that the matching 
conditions are ultra-local in angular coordinates (they do not contain tangential 
derivatives), the most general choice is the following four-parameter family 



<p + = uuacj)^ + uob(f)_ , 

4>' + = uccf)- + u)d(f>'_, (33) 

where 

(f>± = lim <f>(r) , = lim d r (f>(r) . (34) 

r—>R±0 r—fR±0 

Here, u is a complex phase factor, which we include for completeness only. We 
consider real fields and put u = 1. The other parameters obey the restriction 
ad — be = 1. 

There are two important special cases of the conditions (|33|) . Take 

a = d=l, 6 = 0. (35) 

This requires the functions to be continuous across the surface and their deriva- 
tives to have a jump. This is equivalent to having a delta function potential 
V(x) = cS(r — R) on S ^ 1 , which can be viewed as a singular background poten- 
tial concentrated on the surface. The formal limit c — > oo turns this matching 
condition into Dirichlet boundary conditions ((f) ± = 0). 
The other special case is 

a = d=l, c = 0, (36) 

requiring the derivatives to be continuous, and the functions themselves to have 
a jump. This is usually attributed to the presence of a background potential in 
the form of the derivative of the delta function. The formal limit b — » oo turns 
this condition into Robin boundary conditions (| <p± + (f>'± = 0). 

In general, the parameters a,b,c,d,u may depend on the angular coordinates 
on 5 ,£>_1 . In this paper we restrict ourselves to the case where there is not such 
dependence. Then, variables can be easily separated by making the ansatz 

(f) {n) (x) =r 2 -^(f) n j(r)Y {l) (n), (37) 

where Y^(Q) are spherical harmonics depending on the angular coordinates Q. 
Once such ansatz is adopted, the radial functions 4>n,i must satisfy the equation 



d 2 Id v 2 x0 
_j_ _|_ 

dr 2 r dr r 2 71,1 



(31 
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with v — 1 + ^2^, and the matching conditions ( p3[ ) at r = R with shifted values 
of the constants: 

2 - D 1 _ 2 - D , 

a->fl = aH o, c^c = cH a. (39) 

2R ' 2i? v ; 

The degeneracy of each eigenvalue ; is 

v ; /!(£>- 2)! v ; 

In what follows, we will determine the corresponding zeta function and, from 
it, the corresponding heat kernel coefficients. Because we have a continuous 
spectrum, we must separate the translational invariant part (it does not depend 
on the background). We use the procedure described in Ref. |18j using the setup 
of a scattering off the background potential. We have to define the so-called 
regular solutions <ft Pt i which have the same behavior at r — > as the free solution 

<P P A r ) ~ Mpr)- (41) 

The behavior of this regular solution for r — > oo defines the Jost function fi(p): 

Mr) = fi( P )H^(pr) + f;(p)H (1 \pr) . (42) 

In the present case, the eigenfunctions of the Laplace operator can be found 
exactly and they give, for the problem at hand, the Jost function 

fip) = ^[aJ v (pR)HP(pR)' + bpJ v (pR)'HlP(pRy (43) 

- -J u {pR)Hi l \pR) - dJ v (pR)'HV>(pR) 
P 

Now, in order to use a formula like Eq. (§]) we need to have discrete eigen- 
values. So we suppose for a moment that our system is placed inside a sphere of 
larger radius R*. Imposing Dirichlet boundary conditions at r = R*, we obtain 
the following equation for the eigenvalues p = X n /. 

fi(p)H^(pBr) + ftip)H^{pR*) = . (44) 

Then, the (^-function can be represented as a contour integral: 

C(s) = f>0D) f ^-.{p 2 + m 2 r^-\n [ft(p)H®(pir) + fi(p)H<»(pI? 

(45) 

The contour 7 is chosen counterclockwise, enclosing all solutions of Eq. fl44]) on 
the positive real semi-axis and the positive imaginary semi-axis. For convenience, 
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we have introduced an auxiliary mass, which we will later put to zero. There is 
a cut in the complex plane, which goes from im to zoo. Since the number of 
negative modes of the Laplacian is finite, we can always choose m to be suffi- 
ciently large so that all poles of the integrand fl4"5|) are below im. Next, we may 
deform the integration contour as described in |L5], [18|, [T7| to go along the two 
sides of the cut. We perform the limit R* —>■ oo, and drop some contributions 
which are exponentially small in this limit and a term which does not depend on 
the matching conditions (i.e. the "empty space" contribution). The procedure 
sketched above is a quite general one, and not specific of this example, since it 
uses only some general properties of the scattering problem, as Hermiticity and 
ellipticity of the Laplacian. 

Next, we take the limit m — > 0, which is smooth at least for the heat kernel 
asymptotics, and obtain: 

oo 

c(s) = J dk k-^lnf^k)) . (46) 



1=0 



In the Jost function, we can drop any constant factor since it does not contribute 
to (f46D and redefine 



fi(ik) = 1 + jk(IK)' + (3k 2 1' K' + alK (47) 

with new parameters 7 = (d—d)R/ (d+a), (3 = —2bR/ (d+a) and a = 2cR/ (d+a) 
as well as the short-hand notations / = I u (kR), K = K u {kR). 

In order to get the poles of the zeta function (^), thus determining the heat 
kernel coefficients by means of Eq. (|10D, we insert into this Jost function the 
uniform asymptotic expansion, Eq. (0), of the modified Bessel functions and 
obtain 

In /,(<*) = In (l - (3^ d u (t)k d u (t) + a^i v {t)k v {t) + pf(t)k v (t) - i u (t)k d v (t)) 

(48) 

Now, because all functions, i u (t),i*(t), k u (t)k d (t) = 1 + (-) are of order one 
for v — > 00, the leading contribution in the argument of the logarithm is the 
one proportional to (3. As this term grows with u, two cases must be treated 
separately, i.e., (3 = and [3 7^ 0. 

For (3 = we obtain, by means of Eqs. (|55|) , an expansion similar to Eq. 
(|56"D , where the Y^pi are polynomials in the coefficients a and 7. The remaining 
calculations run in the same manner as in the preceding section and we obtain, 
in D = 3 dimensions, 



a x = -4 (a - 7) 7T (49) 
CJ3/2 = (« - 7) 71-2 



13 



2 

a 2 = (5 a 3 - 120 7 - 5 a 2 7 + 3 a 7 2 - 3 7 3 ) 7r 

15 

1 3 

a 5/2 = (a 4 -20a7-2a 3 7 + 367 2 + 2a 2 7 2 - 2«7 3 + 7 4 ) tt5. 
8 

The corresponding results for higher dimensions are given in Appendix 2. 

Next, we turn to the case (3 ^ 0, which corresponds to the presence of a 5' 
potential. Here, we rewrite the logarithm of the Jost function ((471) in the form 



ln/,(zA;) = ln/3 + ln^ (50) 

+ ^ (l + (« - 1) - ^ - - ^ - • 

The first term in the rhs, ln/3, drops out due to the derivative in (|46|) . The 
contributions surviving in the limit (3 —>■ oo in ( |50D are just the same as those 
one obtains for Neumann boundary conditions. Inserting now the asymptotic 
expressions (|55p and proceeding as above, one arrives at the following coefficients 
in D = 3 dimensions, 

Ol = 16— 7T 



13 

„ 3/2 = l(i +2 4| + W ^,,l 



8/3 a 7 1 \ 

a2 = i5U + 6 V +20 ? +8 VJ 1 

a , 7 960 \ 3 . , 

+960 - +m 2. + -y-,. pi) 

Again, the corresponding results for higher dimension are given in Appendix 2. 
As already pointed out in the Introduction, the coefficients present, in this case, 
a dependence on inverse powers of (3. 



5 Conclusions 

In the foregoing sections we used the multiple reflection expansion as a method 
for the calculation of heat kernel coefficient. As a simple example, we considered 
boundary conditions on a sphere and obtained the heat kernel coefficients as a 
finite sum over reflections, Eq. (|32|). An interesting point is a connection between 
Dirichlet and Neumann boundary conditions following from this representation: 
The contributions from the reflections are the same in both cases except for the 
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sign for odd number of reflections. This can be clearly seen already from Eq. (f2~l|) 
and, hence, holds in general. It occurs that this seemingly simple observation had 
not been spelled out before. 

The multiple reflection expansion, as well as the equivalent integral equa- 
tions, Eq. (14) for the propagator and Eq. flTBD for the heat kernel, provide 



a perturbative expansion. For the propagator, this expansion is convergent for 
imaginary frequency (as used in this paper), as was already observed in ||19|| . For 



real frequencies it may diverge. The same holds, presumably, for the heat kernel: 
The corresponding perturbative expansion can be expected to converge. It is 
interesting to note that the convergence of these expansions does not follow from 
a small expansion parameter. For instance, with Dirichlet boundary conditions 
there is no such parameter whereas for matching conditions corresponding to a 
delta function potential on the surface there is one, cf. p2j and Sec. 4. The 



corresponding quantities in the expansion may be numbers which turn out to be 
sufficiently small. In the example with Dirichlet boundary conditions, in Eq. (f28|) 

ujR^- {I v {loR)K u {loR)) < 1 (52) 

holds ensuring the convergence of the geometric series there. 

In general, the convergence issue is not trivial. As an example, we considered 
in Sec. 4 the most general background potential concentrated on a spherical sur- 
face. It is given by the matching conditions in Eq. fl33|), which include a delta 
function potential and its derivative as special cases (Eqs. (|3~5p and (0)). Using 



the techniques introduced in ||15|1 , we calculated for the first time the correspond- 
ing heat kernel coefficients. The lesson with respect to the multiple reflection 
expansion is that, for (3 ^ (in Eq. (|47D (or, equivalently, for b ^ in ( (|33|) )), 
i.e., in the presence of the derivative of the delta function, the coefficients are not 
analytic in (3. In fact, they are polynomials in ^. Hence, the multiple reflection 
expansion cannot converge for (3 ^ 0. 

To summarize, we have stressed the convenience of using the multiple reflec- 
tion expansion for the calculation of heat kernel coefficients, while showing, at 
the same time, some limitations of the method. In general, this method provides, 
after the general Seeley's calculus, the only systematic way to calculate heat ker- 
nel coefficients for manifolds with boundary, and we expect that it will be useful 
in future applications. 
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Appendix 1 



In order to calculate the coefficients an for n > | in Eq. (|32|) it is useful to 
carry out the integration over r in ( |5T| ) using the known formula J dx x I u (x) 2 = 

=j- (^I u (x) 2 ^1 + — /^(x) 2 j. Inserting the second term in the rhs. of Eq. fl23"| ) 
results in the representation 

(53) 

where I = I v {yz)^ K = K v {yz) and we introduced a new variable z = 

Next, we substitute the uniform asymptotic expansions of the modified Bessel 
functions for v — > oo, z fixed: 



s = y - n * 1 ' " ' /_u 



IviyZ) — ^2m^ (l+z 2 ) 1 / 4 *"(^)' K v {yz) — \/~2u (1+z 2 ) 1 / 4 ^(^)> 



(54) 

with 

«t) = *) = E^, (55) 

r>0 r>0 
r>0 r>0 

Here, the notation £ = + z 2 is used. The Debye polynomials u r (t) and f r (£) 

can be found in PBfl , they contain powers of £ from r to 3r. We dont need the 
function rj, since it cancels out in our case. 
We can thus write 

i \ ,M-k 2n 2n+1 p +p+2i-l 

, Hl + j_, _ a Kdr =^_y.j:j: yJ— + - . (56) 

' fc=0 p=fc+l 8=0 

where the coefficients Y^ pi can be calculated easily using a simple computer pro- 
gram. In fact, Eq. fl56|) is the definition of the Y^. For Dirichlet (/i = +1, u — 0) 
and Neumann (// = — 1, tt = 0) boundary conditions they are pure numbers; for 
Robin boundary conditions (/i = -1, m ^ 0) they are polynomials in u. In Eq. 
(|56D the dots denote higher order terms which do not contribute to the considered 
heat kernel coefficients. 

For the integration over z, we use the formula 



r(-s)r(s+ 



2r(^i) v ; 
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One can easily check that the terms with t° from the asymptotic expansions (p5|) 
are cancelled after substitution inside the brackets in K D R (see eqs. ( p8|) and 
(|29|)). This means that any new reflection contributes at least one power of t to 
the integrand in (|53"D. Therefore, according to (|57|) , only several first terms of 
the multiple reflection expansion contribute to any given heat kernel coefficient. 
This explains the finite range of the summations in (|56|). 

The sum over / produces Hurwitz zeta functions. When taking this into 



account, we obtain for the contribution of k reflections to a n (see (|32|)): 



i fi „3/2 2n + 1 P P2s,,fc+1 

s = I - n 1 ^ S J p =fc+ i i=0 z 

The calculation of the residua can be carried out, again, using standard computer 
algebra programs, which leads to the coefficients al in Eq. (|32|). 



Appendix 2 

Here we present the results for the heat kernel coefficients corresponding to the 
matching conditions in Sec. 4, in the cases of some higher dimensional spaces. 
For (3 = we obtain, instead of Eq. ([49]) 
in D = 4 dimensions, 



ai = —2 (a — 7) 7r 2 
1 

2< 



1 1 ^2 4 
«3/2 = -(a -7) 7T 2 



a 2 = -i(2a 3 - 1057 - 2 7 3 ) 7T 2 
6 

a 5 / 2 = — (-3a 2 + 4a 4 -172a7-4a 3 7 + 2987 2 + 3a 2 7 2 - 10a7 3 + 77 4 ) ^ 
64 

in D = 5 dimensions, 

ai = — (a - 7) 7T 2 
a 3 /2 = g (a - 7) 7T 2 



a 2 = (5a 3 - 4657 + 5a 2 7 - «7 2 - 97 3 ) 7r 

45 

a 5/2 = — (-2a 2 + a 4 -76a7 + 1287 2 -4a7 3 + 37 4 ) vri 
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in D = 6 dimensions, 
di = — (a — 7) 7r 3 

1 / s2 I 

«3/2 = -(a -7) 7T2 

a 2 = --^ (6« + 2a 3 - 2977 + 4a 2 7 - 67 3 ) vr 3 



1 7 
— (a 2 + 4a 4 - 508«7 + 4a 3 7 + 802 7 2 -a 2 7 2 - 26«7 3 + 197 4 ) Tr* 



a 5/2 = 12g 
and in Z) = 7 dimensions, 

ax = -— [a -7) vr 
15 

4 2 1 

^3/2 = — (" -7) TT 2 



a 2 = __ (5« 3 - 10507 + 15a 2 7 + 3«7 2 - 237 3 ) vr 



1 7 
a 5/2 = — (-6a 2 + a 4 -172a7 + 2a 3 7 + 2807 2 - 10a7 3 + 77 4 ) 7T2. 
30 



For (3 7^ we obtain, instead of Eq. (pip, 
in D = 4 dimensions, 

8 2 
01 = r 



If a 7 128\ 5 

«s/2 = -^9 + 64--3 2 l + -)^ 



16 / a 4 , 

a2 = YlV + ^ )7r 



2 



^(-59-512^-2241 -512^ + 4096^ + 2048^ + 5122; 

5 



Ot 1 1 

+32768 — + 8192 + 32768— ] ~" 



in D = 5 dimensions, 

32 vr 2 



«1 



3 /3 



2 / a 7 4 \ 5 

a 3/2 = - 3 + 2 2 - + — 7T2 

3/ 3 V P P P 2 J 

16/1 a 7 80 . 

° 2 = 45^^ + 6 V -2 V + ^ )7r 



2 



4/17 a 2 a 2 a 16 
3 V 240 ~ p ~~ p 2 + p 2 + ]P + p* 



5 

7T2 
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in D = 6 dimensions, 



4 3 

1 /191 , a 7 128\ 

2/3 a 7 16\ o 
02 = sU^V"^^ 



^ = 12188 ( 103 + 512 I- 48 °l + f + 12288^- 6144 |I "512 J 



a 7 98304\ 7 

+98304^-24576^ + — W> 



and in D = 7 dimensions, 



3 



64 7T 

01 = 15^ 

/ a 7 16 . 

a 3/2 = 4(l0 + 8--16-^ + -)7r5 



4/15 „ a „ 7 12 a 2 cry 

as/2 = islr^-VV^V F 

a 7 32\ 7 
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